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Abstract 

A large amount of population models use the concept of a carrying capacity. Simulated populations are bounded 
by invoking finite resources through a survival probability, commonly referred to as the Verhulst factor The fact, 
however, that resources are not easily accounted for in actual biological systems makes the carrying capacity parameter 
ill-defined. Henceforth, we deem it essential to consider cases for which the parameter is unnecessary. This work 
demonstrates the possibility of Verhulst-free steady states using the Penna aging model, with one semelparous birth 
per adult. Stable populations are obtained by setting a mutation threshold that is higher than the reproduction age. 
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1. Introduction 

The evolution of a population can be described using 
the difference equation 



M+i = rN,\ 1 - -^ 



(1) 



where A^, is the population at time t, r is the intrin- 
sic relative growth rate, and K is the carrying capacity. 
This logistic type of growth was first introduced by Ver- 
hulst in the mid 1800s 1 1]. In an environment where re- 
sources are finite, populations are bounded. The equiv- 
alent differential equation leads to a final state which 
is a fixed point. Discrete counterparts, such as the one 
presented above, result in bifurcations, limit cycles and 
chaos |l, 2|1. An in-depth discussion of the more com- 
plex set of solutions obtained by considering various as- 
pects of the Verhulst equation may be found in [3]. 

The concept of a carrying capacity is commonly in- 
corporated in population dynamics models. Limitations 
on food supply, space and other necessities suppress 
growth. Organisms compete in order to stay alive. Each 
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subject to a survival probability, V, = 1 - N,IK, also 
known as the Verhulst factor ||4|, |5|] . The main moti- 
vation for its introduction is to make certain that pop- 
ulation sizes are kept finite \^2t^Sl- This parameter 
may be applied as a constraint for simulated populations 
in order to maintain their size within computer limits. 
However, the carrying capacity is not constant llOll . An 
environment's actual ability to support life highly de- 
pends on prevailing knowledge and technologies fll|]. 
It is not easily determined 1,12] . neither is it well de- 
fined. 

We choose to consider cases for which the carrying 
capacity is not essential. Because resources are not ex- 
plicitly accounted for in many observed systems, we be- 
lieve this case can be of great relevance to biological 
population studies. 



2. The model 

The Penna model fl3J] is apopular technique for sim- 
ulating aging populations \\% . Individual characteris- 
tics are represented by a string of binary numbers of 
length L. Over time, mutations accumulate thereby 
modifying an organism's viability, survival and fertility 
OlSll . For each year added to an individual's age, one bit 
in the genome is read. Active genetic traits, therefore, 
are located at bit numbers (locus) less than or equal to 
the current age. Zeroes correspond to healthy genes and 
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ones are bad or mutated genes. An individual dies a ge- 
netic death when it reaches old age (age - L) or when 
the number of expressed harmful mutations equals the 
threshold value, T. Here, only bad mutations are con- 
sidered because helpful ones, in comparison, are very 
rare lll3ll . Reproductive maturity is achieved at age R. 
Starting from this age until death, individuals generate 
B offsprings per time interval, where B is the birth rate. 
In asexually reproducing populations, newborns copy 
the genes of the parent (both active and inactive) and 
acquire M additional mutations. The mutation rate, M, 
can take on any value but is usually set to one since ex- 
pression of new mutations is not very frequent in nature 

To include the effects of finite resources, individuals 
are subject to the same survival probability regardless 
of the quality of genome, fitness and age - VA imple- 
mentation. It has been argued [17.1 that older individ- 
uals should be more adaptable and, thus, have greater 
chances of survival. Furthermore, it was added that due 
to the Gompertzian behavior of the original model, lim- 
iting the number of births is enough to avoid exponential 
growth. Hence, an alternative approach was proposed 
applying the Verhulst factor to newborns only (the VB 
implementation). For lack of a vital role in popula- 
tion dynamics |6], efforts have been made to eliminate 
the random deaths imposed by the Verhulst factor. In 
ILZI], a Verhulst- free Penna model was implemented on 
a square lattice, where each lattice site can accommo- 
date only one individual at a time. The abstract limit 
on the capacity of the environment is, thus, given by the 
size of the grid. Another study [18j modified the model 
such that for every death, an adult is randomly chosen to 
produce one offspring (a replacement). This procedure 
keeps the total population size constant. 

In this work, we modify the 32-bit Penna model, 
making it Verhulst-free. Deaths are now due only to 
genetics - old age and mutation accumulation. To limit 
growth, we allow individuals to breed only once in their 
lifetime (at age R). This behavior is observed in semel- 
parous populations such as the mayflies and the Pa- 
cific salmon lll9l l20l 12111 . The simulation starts with 
No = 20000 perfect individuals (no bad genes) and runs 
for 3000 iterations. The choice of evaluating deaths be- 
fore births limits the maximum reproduction age value. 





Figure 1: Exponential decay curves, NIN„,ax = e "', associated with 
T < R: (A)T = l,aaye = 0.081, (B) T = 5,0,,,,, = 0.015, (C) 
T = 8,aa„ = 0.011 and (D) T = W,a„„c = 0.006. Extrapolation 
shows that extinction begins at later times for populations with higher 
reproduction ages, R. 



parameter and the reproduction age. We observe expo- 
nential decrease when T < R. Stability, on the other 
hand, is achieved by setting T > R. The populations 
presented hereafter are normalized about the maximum 
value obtained for each parameter set. 

3.1. Extinction cases 

Extinction has two causes - mutational meltdown and 
overpopulation [22]. The latter results from a sudden 
fluctuation in population size beyond the carrying ca- 
pacity. But without the Verhulst factor (no stochas- 
ticity), the only possible causes of death are old age 
and mutation accumulation. At low T, selection favors 
deaths at younger ages in order to keep the gene pool 
clean. Thus when T < R, a good fraction of the pop- 
ulation dies prior to reproductive maturity. With births 
limited to one offspring per individual, the number of 
newborns is not enough to compensate for losses due to 
bad mutations. The population dies of mutational melt- 
down. 

Best-fits show delays in the onset of extinction for 
those that reproduce late in life, large R (Fig. [TJ. The 
rate of decay is effectively slowed down by an increased 
tolerance for bad mutations. Regression analysis yields 
a power law relationship between the mutation thresh- 
old and the average decay rate (Fig.O. 



3. Analysis of results 

In the B - I limit, the behavior of the resulting popu- 
lation is dependent on the relative value of the threshold 



3.2. Nonzero steady states 

Stable populations were verified using first return 
maps 12311 . Consider the sequence jcq, xi, X2, X3 . . . x,_i 
and Xf. The first return map is a plot of consecutive 
entries of this sequence, i.e., (xq, xi), (xi,X2), (x2,xj,), 
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Figure 2: Average decay rate associated witli R = 14,21 and 28. The 
rate of extinction decreases with increasing tolerance for harmful mu- 
tations. The relationship is a power law, y = 0.08.V"' . 



. . .(x,_i,x,). Fig. [3] presents the return maps of the 
steady state populations associated with R - 4,^ and 
12 when B = 1 and T - 15. The number of fixed points 
reflect periodic variations in the time series. For these 
steady states, we find that the magnitude of the fluctua- 
tions in total population size increases with R. 

To minimize fluctuations, running averages were ob- 
tained for the time series. This was done using inter- 
vals of 50. Variations in the transient part of the curves 
(FiglH) indicate differences in saturation time. The slope 
from the highest value to the onset of steady state be- 
comes steeper as R increases. Hence, those with higher 
reproduction ages take longer to equilibriate. Similar 
to the original Penna model results [13], steady state is 
achieved faster by allowing more mutations, that is, by 
increasing M. 

Fig. |5] shows the variation of the mean population 
with threshold, for a given R. This corresponds to the 
average population size of the last 300 iterations of a 
particular simulation run. At these times, the simu- 
lated populations have already achieved steady state. 
Note that the mutation threshold can have values from 
1 to L, the bit-string length. However, since stability is 
achieved only when T > R (for the B - 1 case), the 
number of T values that result to nonzero steady states 
decreases as R increases. We find that the total popu- 
lation scales with threshold. Considering all possible R 
values (from 1 to L - 1), we obtain linear fits that have 
slopes ranging from 0.015 to 0.035. 

The steady state associated with T = L (or M = 0) 
takes advantage of the finite bit-string characteristic of 
the Penna model. Following the derivation in fl9], we 
let A^(fl, t) be the number of individuals with age a at 
time t. The limit on the length of the string forces deaths 




Figure 3: First return maps of stable populations resulting from B = 
1, r = 15: (A) R = 4, (B) R = 8 and (C) R = 12. These correspond to 
the last 300 iterations of the simulation, at which time the populations 
are already at steady state. 
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Figure 4: Normalized populations associated with B = 1,T = W: (A) 
R = 3, (B) R = 6 and (C) R = 9. Steady state is achieved faster for 
larger values of the reproduction age, R, and mutation rate, M. 



at age L. In the absence of the other death factors, the 
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Figure 5: Mean normalized populations (circles) obtained from differ- 
ent threshold values and reproduction ages. Linear fits (dashed lines) 
yield (A) y = 0.029a- + 0.052 for R=i, (B) y = 0.031a- for R = 8, (C) 
y = 0.031A + 0.009 forR = 12 and (D).y = 0.019a + 0.076 forR = 24. 



discrete time evolution of the population is described by 
N{a, t) I <a < L-l 



Nia+l,t+l) 







fl > L- 1 



and 



N{0,t+l)^BN(R,t+ 1). 



(2) 



(3) 



At steady state, the average number of newborns is con- 
stant, i.e., A^(0, t) - N{Q). From the equations above, we 
have 



N(R, f -H 1) = N(0, t-R + l)^ N(0). 



(4) 



Using (HJi in ^, it is straightforward to show that in 
the no mutation case (T - L or M = 0), stability is 
achieved only when B - I. No other nonzero steady 
states were found. Beyond B - I, populations either 
increase indefinitely or tend to zero (Fig|6]l. 

It has been shown ll24ll that asexual populations have 
at most r - 1 mutations before the minimum reproduc- 
tive age. For semelparous species, all bits above age R 
are set to one. In [25], we have demonstrated that this 
is not always the case. As T increases, the strength of 
selection decreases and more mutations are preserved in 
the genome. When T > R, all bits at steady state contain 
defective genes. The process of mutation accumulation 
is very fast in asexual populations. To slow down the 
process, one may consider introducing recombination, 
as in sexual populations, or allowing positive mutations. 

4. Summary and conclusion 

It is common practice in population models to at- 
tribute stability with the concept of a carrying capacity. 



Figure 6: Population explosion curve associated with the Verhulst- 
free Penna model. Under the curve, behavior other than exponential 
decay is observed only if B = 1 and T > R (nonzero steady state). 
All threshold values not included in the figure result in exponentially 
increasing populations whenever B > 1 . 



Verhulst-free cases generally lead to Malthus catastro- 
phe - exponential growth or decay. The carrying ca- 
pacity, however, is an ill-defined parameter 1 10, 11, 12|]. 
It is for this reason that we explore cases for which the 
concept is unnecessary, that is, to find stable populations 
without having to impose a Verhulst factor 

Within the framework of the Penna model, the artifi- 
cial cap on the individual lifespan imposed by the bit- 
string length and the semelparous limit on reproduction 
make possible the Verhulst-free implementation. 

Population explosion and extinction are still observed 
in many of the Verhulst-free cases. We find those asso- 
ciated with B - I and T > R result in nonzero steady 
states. 
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